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We study the static potential of a color singlet quark-antiquark pair with (fixed) distance r in D = 3 
and D — 2 space-time dimensions at weak coupling (ar C 1 and gr <C 1, respectively). Using the 
effective theory pNRQCD we determine the ultrasoft contributions, which cannot be computed in 
conventional perturbative QCD. We show in detail how the ultrasoft renormalization in pNRQCD 
is carried out. In three dimensions the precision of our results reaches 0(a 3 r 2 ), i.e. NNLO in the 
multipole expansion, and NNLL in a a/AV expansion, where AV ~ aln(ar). We even present 
results up to partly N 4 LL order and compare them to existing lattice data. Finally we discuss the 
relevance of the perturbative calculation in two dimensions, where the exact result is known. 

o ■ 

CN ' PACS numbers: ll.30.Pb, ll.10.Hi, 12.39.Hg 

CD 

fj^ ■ I. INTRODUCTION 

OS 

The singlet potential between two color sources in the static limit, i.e. with infinite masses, is one of the 
most accurately studied objects in QCD. It plays an important role in understanding the dynamics of QCD. 
On the one hand it is an essential ingredient in a Schrodinger-like description of heavy quarkonia. On the 
other hand, a linear growing potential at long distances (in the large N c limit) indicates confinement of the 
strong interactions. This has motivated a lot of effort to steadily increase the accuracy of the perturbative 
prediction for that object [H-[l0j. which should be reliable at short distances. A crucial point is obviously 
to distinguish which region can be described within perturbation theory and to understand the crossover to 
the non-perturbative regime 

una, a- 

Besides the singlet also the color octet static potential has attracted some interest in the last decade. It has 
I/-} been computed with two-loop precision at weak coupling in Ref. UM . Together with the so-called gluelumps, 

it describes the short-distance behavior of the hybrid potentials [15[ [l6[ , which are potentially important for 
the theoretical prediction of physical hybrids made of heavy quarks. At long distances the hybrid potentials 
could also help to shed light on the question whether the dynamics responsible for confinement is of string 
{SJ ■ type and if so, to determine its structure exactly. 

In this paper, we will focus on the singlet static potential, V 8 (r), and on the singlet static energy, E s (r), 
at short distances. V s (r) is the leading order potential introduced in a Schrodinger-like description of 
heavy quarkonia. E s (r) is the energy of a static color singlet quark-antiquark pair with (fixed) distance r, 
being the object actually computed on the lattice. Its physics at short distances is governed by at least 
two physical scales. One is the soft scale ~ 1/r, and the other one is the ultrasoft scale ~ V s . The most 
convenient framework to perform the corresponding calculations is the effective field theory (EFT) "potential 
NRQCD" (pNRQCD) [lj] (for a review see [Hj]). pNRQCD also has been instrumental in many of the latest 
perturbative results quoted above, especially in the determination of the ultrasoft contributions. In its static 
limit the theory exploits the scale separation 1/r ^> V s . V s accounts only for the soft scale contributions 
to the static energy and can be computed in perturbative QCD (pQCD). E s also includes effects at the 
ultrasoft scale, which in the static limit is generated by a resummation of certain classes of loop diagrams to 
all orders and serves as an infrared (IR) regulator [19| . These contributions can be calculated in pNRQCD 
from a finite number of diagrams. 

In order to gain a deeper understanding of the static potential one can consider how it is qualitatively 
affected by changing the number of dimensions from four (4D) to three (3D) or two (2D). The three- 
dimensional result is also important on its own. Four dimensional thermal QCD undergoes effectively 
a dimensional reduction for large temperatures. Therefore, determining the renormalization group (RG) 
structure for the static potential in three dimensions may open the way to a resummation of logarithms 
at finite temperature. Three dimensional space-time is moreover a good testing ground for renormalon 
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issues, since the linear power divergences associated with renormalons in four dimensions become logarithmic 
divergences in three dimensions and can be traced back using dimensional regularization. The computations 
in three and two dimensions finally represent consistency checks of the theoretical approach used to describe 
the 4D potential. 

Partially driven by these motivations the computation of the static potential has been carried out up to 
two loops for arbitrary dimension D in Ref. 0|. Whereas in four dimensions the result is finite, this is not 
the case in three dimensions. Beyond one-loop it suffers from IR singularities, whose origin was obscure 
at the time they were encountered. In particular, it was not clear whether the space-like end-strings of 
the rectangular Wilson loop, which defines the static potential non-perturbatively, could contribute, raising 
doubts about the independence of the result on the precise form of the strings. 

In 4D such IR singularities first appear at three-loop level and were shown to cancel with the ultrasoft 
ultraviolet (UV) singularities in Ref. [3|. Up to now an analogous ultrasoft computation has been missing 
in three dimensions. Therefore, the precision of the 3D static potential has been limited to 0(a 2 r) from the 
finite tree-level and one-loop results. With this work we will show for the color singlet state that the UV 
divergences of the ultrasoft pNRQCD calculation in fact cancel the soft pQCD IR singularities exactly. This 
allows us, for the first time, to go beyond the one-loop level and to obtain the logarithmic corrections at two 
loops as well as the full analytical structure of the 0(r 2 ) correction. We will point out the peculiarities of 
the 3D computation in detail. 

The outline of the paper is as follows. In Sec. |TI]we introduce the theoretical setup for our calculations, 
we define the power counting of the EFT in D < 4 dimensions and review existing results for the relevant 
ultrasoft loop diagrams. Sec. Mil contains the complete RG improved ultrasoft calculation up to NNLL and 
partly up to N 4 LL order in three dimensions as well as a comparison of the results to lattice data. In Sec. IIVI 
we discuss the situation in two dimensions and in Sec. [V]we present our conclusions. 

II. PNRQCD 

Up to the next-to-leading order (NLO) in the multipole expansion (and irrespectively of the space-time 
dimension) the effective Lagrangian density of pNRQCD in the static limit takes the form [llj , [HI : 



C us = Trjst (id ~ V s (r)) S + O* {iD - V a {r)) O j 



+ 5 VSi(r)Tr{O t r-ES + S t r-EO} + .g-^^Tr{O t {r -E,0}} +C(r 2 ). (1) 

We define color singlet and octet fields for the quark-antiquark system by S = S(r, R, t) and O a = O a (r, R, t) , 
respectively. R = (xi +X2)/2 is the center position of the system. In order for S and O a to have the proper 
free-field normalization in color space they are related to the fields in Eq. ([1]) as follows: 

1 T a 

C -.S, = ^=O a , (2) 



where Tp = 1/2 for the fundamental representation of SU(iV c ). All gluon and scalar fields in Eq. ([TJ are 
evaluated in R and the time t, in particular the chromoelectric field E = E(R, t) and the ultrasoft covariant 
derivative iD O = id$0 — g[A (Il, t), O]. V/ a ,o,A,B} ( r ) are the Wilson coefficients of the effective Lagrangian. 
They are fixed at a scale v smaller than (or similar to) 1/r and larger than the ultrasoft and any other scale 
in the problem by matching the effective and the underlying theory, which in this case is QCD in the static 
limit. 



A. Power counting 

Because the mass dimension of the coupling is [g 2 ] = M A ~ D , we have (at least) the following physical 
scales involved in the problem: 

k ~ 1/r (soft) , k ~ V (ultrasoft) , k ~ g 4_D (non-perturbative) . (3) 
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In order for perturbation to make sense at the (soft) matching scale we demand 

gr^ «1, (4) 

i.e. weak coupling. Therefore g 2 r A ~ D plays the role of a dimensionless expansion parameter. It is a good 
expansion parameter at short distances for D < 4 (for D = 4 Eq. (U) implies g <C 1). 

We are left with two scales: V ~ 1/r x g 2 r 4 ~ D and g^-° much smaller than 1/r as we go to short 
distances. In principle, this guarantees that the multipole expansion makes sense, but it says nothing about 

2 

the relative size of V and g 4 -° , i.e. whether one can use perturbation theory at the ultrasoft scale V. On 
the other hand note that 

V / g ^o „ r 3-D g 2(i^) „ ( rg ^o f-D (5) 

Therefore, D = 3 is the turning point for the use of perturbation theory at the ultrasoft scale. 
For 4 > D > 3 we have a perturbative expansion in g: 

1/r > V > g*=v . (6) 
In this case there is a double expansion with the two parameters 

2 

For _D < 3 the right expansion parameter would be V/g 4 - D . In principle, this is not a real problem. 
We just have to expand the exponent in the ultrasoft correlators, as we will see later, while the multipole 
expansion is still valid. This is potentially quite interesting. In two dimensions the exact result is known 
and equivalent to the tree-level result. Therefore, strong cancellations are supposed to occur among the 
contributions from the different scales: soft, ultrasoft and g. We will further discuss this issue in Sec. IIV1 

The case D = 3, the main focus of this paper, deserves a special discussion. At short distances we find 
V ~ g 2 \n(rv). This implies that 

g*lV.~\/Up,) <1, (8) 

if v ~ V and we conclude that we can use perturbation theory at the ultrasoft scale V. We therefore formally 
distinguish between the scale V and g 2 . Logarithms from the ultrasoft perturbative computation will have 
the form IniV/v) and are rendered small, if we set v ~ V. Thus it is legitimate to consider the ultrasoft 
regime as perturbative, i.e. the pNRQCD loop expansion makes sense (for sufficiently small r). The inclusion 
of non-perturbative effects will be discussed in subsection IIIIDI 



B. Bare data 



Here we summarize some bare results that will be relevant in the next sections. In the following we will 
use the index U B" to explicitly denote bare quantities. Parameters without this index are understood to be 
renormalized. 

The general expression for the bare singlet potential in momentum space, V s ,b, ™- D dimensions can be 
written as 

^-^I^jl + E^k-^M^j . (9) 

The coefficients C\{D) and 02(D) can be found in Ref. |2|. Throughout this paper we will use the notation 
D = 4 + 2e = n + 2e„, where e„ = D 7, n parametrizes the (typically infinitesimal) difference to the closest 
integer dimension n = 4, 3, 2. 

After the Fourier transformation to position space Eq. ([9]) becomes (see for instance Ref. [Ic| ) 



V S ,B - -C F g B > 2-2» f ,3/2+^n y W > 



2 V- 5l"^ 2( " +1)e Cn{D) r[l/2 + {n + l)e] 
1b 

-Cf9% 



n=0 



(4 7r )nD/2 22-2n6 7r 3/2+er[l 



g 2 B n c n (D)r- 2 ( n +^ 



n=0 
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The singlet static energy can be considered to be an observable for our purposes. It consists of the 
potential, which is a Wilson coefficient, and an ultrasoft contribution 1 , either bare or renormalized: 

E s (r) = V S , B + 6E™ B = V sm + SE™ m . (11) 

The ultrasoft contribution can be expressed in a compact form at NLO in the multipole expansion (but exact 
to any order in a) through the chromoelectric correlator. It reads (in the Euclidean) 

rp />oo 

SKb = V A {D J l)N f I dte- 1(% - B -^ fl) <«ac| 5 ES,(*)^(t 1 0) ff E|,(0)|i;ac) . (12) 

Concrete results for the ultrasoft corrections in D dimensions are known at one loop since Ref. [2l| (see 
also 0,i]). The pNRQCD one-loop computation yields 

6E™ B (1 - loop) = -g 2 C F V 2 (l + £) r[2 + e]T[-3-2e] rA ^ +2£r ^ 

where we have defined AT^ = V — V s . The two- loop bare expression can be deduced from the results obtained 
in Refs. [3, [22| and reads 

5E? B (2 - loop) = g 4 C F C A V 2 T[-3 - 4c] \v^(e) - (1 + 2e)v[ 1 \e)] r AFj +4e r , (14) 



where 



p(1) ^)= (2 ^ )2 47r2 1 +2£ r 2 [i + £ ] g ( e ), (15) 



^ 1)(e)= (2^4 7 r^ r2[1 + £]gl(e) ' (16) 



and 



2e 3 + 6e 2 + 8e + 3 _ 2 e r(-2e - 2)r(-2e - 1) 
9[e) ~ e(2e 2 + 5e + 3) (2e + 3)r(~4e - 3) ' ( ' 

_ 6e 3 + 17e 2 + 18e + 6 4(e + l)n f T f 2 (e 2 + e + l) r(-2e - 2)r(-2e - 1) 
31 W ~ e 2 (2e 2 + 5e + 3) + e(2e + 3)N C + e(2e + 3)r(-4e - 3) ' ( ' 

C. Renormalization, generalities 

The bare parameters of the theory are as (<7b) and the potentials Vb = V{ s ,o,a,b},b- I n our convention 
as has integer mass dimension and is related to gs by 

n 2 u 2t " 

a - = £ ^-» (19) 

where v is the renormalization scale. It has a special status since it does not receive corrections from other 
Wilson coefficients of the effective theory. Therefore, it can be renormalized multiplicatively: 

a B = Z a a , (20) 



1 If one has enough precision non-perturbative effects at the hadronization scale Aqcd ~ 9 4 ~ D should also be included. We 
will discuss them in Sec. I1I1I 
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where 

oo 

Z^l + ^Z^-. (21) 



8=1 r ' 



The renormalization group equation (RGE) of a is 



In the limit e„ — > 



i/— a = a/3(a; e n ) = 2e„a + a/3(a; 0) . (22) 



v^-a = ap(a; 0) = a /3(a) = -2a— Z^ . (23) 
dv da 



The bare potentials Vb in position space have integer mass dimensions (note that this is not true in 
momentum space) and, due to the structure of the theory, we do not renormalize them multiplicatively, see 
the discussion in Ref. M- We define 

V B = V + SV . (24) 

5V will generally depend on the (matching) coefficients of the effective theory, i.e. on a and V and on the 
number of space-time dimensions. In D(n) dimensions, using the MS renormalization scheme, we define 



oo 1 

5V = Y J Z { y ) -. (25) 



s=l 6,1 

From the scale independence of the bare potentials 

v JLv b =0, (26) 
dv 

one obtains the RGE's of the different renormalized potentials. They can schematically be written as one 
(vector-like) equation including all potentials: 

V ±-V = B{V), (27) 



B(V) = - {u^w) . (28) 

Note that Eq. (|27l) implies that all the l/e„ poles disappear once the derivative with respect to the renor- 
malization scale is performed. This imposes some constraints on 8V: 

0(l/e n ) : B(V) = -2a— , (29) 

oa 



t'i 1/4) : fl(F) A Z W + ap(a)-^-Z^ + 2a-^zf = , (30) 



and so on. 



III. PNRQCD (D=3) 

In three dimensions the purely gluonic sector is superrenormalizable. As a consequence the coupling 
constant is not renormalized: 



ub = a 



(31) 



G 



=0= ® = 
FIG. 1: One-loop contribution to the octet propagator. The dotted line represents the A field. 

At leading order in the multipole expansion the singlet field of the quark-antiquark system is free, i.e. it 
does not interact with gluons, and renormalization scale independent. In other words its renormalization 
constant (S_b = z]^ 2 S) is 

Z s = l + 0(r 2 ). (32) 
Similarly the singlet potential is not renormalized at 0(r ): 

SV S = 0(r 2 ) . (33) 

For the octet field the situation is different. Even at leading order in the multipole expansion it has 
a residual interaction with ultrasoft gluons. The octet potential receives an ultraviolet (UV) divergent 
correction from the one-loop self-energy diagram shown in Fig. [TJ It gives 

Higher loop diagrams cannot contribute at O(r ). This is because the coupling and/or the potentials 
must appear perturbatively (with positive powers) in the Z's. Since a has positive mass dimension, the 
potentials would appear with negative powers in higher loop corrections to Eq. (|34p , which is not allowed by 
renormalizability. By the very same reason the octet field does not require renormalization at 0(r ), i.e. 

Z = l + 0(r 2 ). (35) 

This will be all we need to know about the renormalization of the octet sector. 

The operators that appear at 0(r) in the Lagrangian, Va/Bi are not renormalized at 0(r ) either, i.e. 

Z A/B = 1 + 0(r) . (36) 

The reason is, as mentioned before, that the coupling and the potentials have to appear perturbatively in 
the counterterms. Otherwise the renormalizability of the theory at LO of the multipole expansion would be 
spoiled. Moreover, we can take Va — Vb — 1, as 0(a) soft corrections would be multiplied by factors of r 
and would move us away from the precision of V s aimed for at this paper. Actually, from inspection of the 
possible diagrams that will contribute at the soft scale, we know that Va/b = 1 + ^(a 2 ) |7|. 

We now focus on the singlet and the renormalization of V s beyond O(r ). The singlet potential is IR safe 
up to soft one-loop order. At two soft loops in dimensional regularization IR poles up to 0{l/e\) appear Q. 
The ultrasoft computation in pNRQCD yields the following results for the counterterms 



'^-r 2 AV x C A a 2 C F V 2 ^ 
(13tt 2 - 2208) C\ + 8 (19tt 2 + 144) C A T F n f - A8T F n f (4 (tt 2 - 10) C F + Tr 2 T F n f ) 



Z$ iX = -r 2 AV 2 aC F V 2 --r 2 AV x C A a 2 C F V 2 - (37) 



2304 



zf a X = ~r 2 AV x a 2 C F C A V 2 - - r 2 a 3 C F C 2 A V^ , (38) 



Zg> =-r 2 o?C F C 2 A V 2 A ^ (39) 
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Z ( ™J = V n > 3 . (40) 

These are the complete C(r 2 ) results. We used the MS renormalization scheme to derive them and 
we checked that their form is invariant under scheme transformations that amount to a change of the 
renormalization scale. In other words the form of Zy o and Zy a remains unchanged, if we replace a — » qcj 3 , 
where cx is a 63-independent constant (cms — !)• We refer to this class of renormalization schemes as 
"global" and indicate globally renormalized quantities by an index "X" representing the scheme (X=MS, 
MS, etc.) as e.g. AVx in the above equations. 

The fact that one can renormalize the potential with a finite number of terms at a given order in the 
multipole expansion (at 0(r 2 ): Zy = for n > 3) reflects the super-renormalizability of the theory. Let us 
now explain how Eqs. Q37H39P are obtained. 

1) The first term in Eq. (|37|) comes from the 1 /e$ divergence of the ultrasoft one-loop correction in Eq. (|13[) . 
It has a scheme independent form and fixes, together with Eqs. (j2U)l and (|3U|) . the first term of Eq. (|3"5)l 
and Eq. (|3U)) . Of course it is also possible to compute these 1 /e§ and l/e\ divergences directly. They are 
generated by the entanglement of the ultrasoft one-loop diagram with the octet self energy in Fig. [Tj 2 Other 
types of diagrams with the same number of loops are less divergent. 

2) The second term in Eq. (|37p follows from the remaining l/e3 divergence in the ultrasoft two-loop 
computation, Eq. (fl"4"|) . once all subdivergences (associated with the octet potential) have been subtracted. 
We have checked its (global) scheme independence explicitly. This result combined with Eq. (|3"0")) then fixes 
the second term of Eq. (|38)l . 

3) We now compare the results of 1) and 2) (with Va = 1) to the two-loop calculation of the soft 
contribution, which equals the bare static potential of Ref. ■ Keep in mind that the comparison is carried 
out in D = 3 + 2e3 dimensions and AV is a polynomial in €3. We find perfect agreement of the soft IR 1 /e§ 
and 1 /e§ singularities and the corresponding ultrasoft counterterms in Eq. (|38| and (|39|) . 3 That means that 
the respective divergent parts of the bare quantities in Eq. ([Til cancel. Note that this can be understood as 
a non-trivial check of two independent determinations of these terms. 

4) Our findings so far open the way to obtain the missing subleading UV divergence in Eq. (f37|) . The 
key point is that the counterterms of the potential cannot contain negative powers of AV. They would give 
rise to negative powers of \n(rv), which cannot be absorbed by the potential (cf. Eq. (ITU1) '). Hence Zy^ 
is a polynomial of AV. From the previous computations we know all terms of the polynomial except the 
constant at O(AV ). This term can be inferred from the (soft) result in Ref. [2J. We simply subtract all 
the l/e| and l/e§ counterterms as well as the first two l/e^ terms according to Eq. (|37l) . The remaining 
(IR) 1/63 divergence must then match the last term in Eq. (f3"T|) in order to render Eq. ([TTT) finite. In this 
way we fix the AV° piece of Zy indirectly. A consistency check of this computation is that the result is 
really independent of AV (i.e. free of ln(r)) and invariant within the class of global renormalization schemes 
defined above. Thus we have obtained the complete renormalization structure of V s at 0(r 2 ). 



A. pNRQCD RG 

We will now derive the RG evolution of the matching coefficients of the effective theory. The running of 
the coupling constant is trivial, since it is not renormalized in three dimensions, as follows from Eq. (|31l) : 

V—a = + O(e 3 ). (41) 

As a consequence of Eq. (f3T>l) Va/b docs not run either at O(r ). Similarly there is no running of the singlet 
static potential at LO in the multipole expansion, whereas the running of the octet potential according to 
Eqs. (IMl) and ([34]) reads 

v^-V = -C A a. (42) 



2 In fact one can even perform an exact resummation of these diagrams (at least of the divergent pieces) to handle all those 
divergences at once. 

3 Actually, if we compare with Eq. (3.51) of Ref. we find a discrepancy, but if we explicitly sum up the expressions for the 
contributing diagrams given in the same reference the results agree. Most likely there is a typo in Eq. (3.51) of Ref. Q. 
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We are interested in the running of the octet potential only insofar as it enters the running of A.V, which 
including the tree-level matching condition is then given by 



AV x (r; v) = V , x (r; v) - V s>x (r; v) = -aC A \n(rud x ) + 0(r 2 ) + 0(e 3 ) . 
The index X stands for the scheme (e.g. MS or MS): 



(43) 



d MS = e 1E/2 ^ ~ 2.36546, 



''MS 



d-MS < 



MS 



e~> E /2 ~ 0.890536, 



(44) 



with 



MS 



e i/2(in(4ir)-7 E ) This result will be sufficient for our purpose, which is the computation of the 
singlet static potential and energy. 

From the counterterms determined in the previous subsection, we can derive the complete running of the 
singlet static potential at 0(r 2 ): 



V ^V., X = ^fvy(AV x ) 2 a + C F Vir 2 AV x C A a 2 



(45) 



r 2 a 3 C F Vl 



(13tt 2 - 2208)C^ + 8(19tt 2 + lU)C A n f T F - 48n / T F (4( 7 r 2 - 10)C F + Tr 2 n f T F ) 

384 



Its form is globally scheme independent. 

By solving the RG equations we obtain for the 0(r 2 ) contribution 



where 



V s M *(v) = V s MS (r;v = -) + K«^ MB (r; v) , 
r 



(46) 



V s RG -> MS (r;v) =C F r 2 a 3 
' 1 



x -Ciln J M 



1 



-C 2 A [-2 + "f E + lmr] \n 2 {rv) 

c ' 2a ( w + AlE ~ 46 + {2lE ~ 4) ln7r + ln2 



+n f T F (c A (3 + 2j£) + C F (5 - ) - (n f T F ) 2 £ 



ln(rz^) 



(47) 



is the running and 



T/ MS (r;i/=-) = C F a\n{r 2 v 2 s ^ E ) + -C F (7C A - 4n f T F )a 2 r 



' tt 2 , 209 43 , , 

^(39 7is - 715 + 564 ln 2 + 39 ln^) + — - -£(3) 

1e - q 1e - 138 7b + 3 (7b - 47b - 46) ln tt + 3(7 B - 2) In 2 vr + In 3 tt 



4.s 



r 7r 1 
-n f T F C F —(15 - 6 7B - 81n2 - 61htt) + ~(5j E - 8 + 51n7r) 



-n f T F C A 
-(n f T F ) 2 



.24 

TT 2 1 

— (57 7 _e - 97+ 121n2 + 571n7r) + -(97b - 22 + 9 ln tt + 10C(3)) 
288 6 



48 



(7 - 3 7B - ln(16^ 3 )) 



(48) 



is the initial matching condition. Note that the tree-level and one-loop matching conditions have been 
included. For the two- loop initial condition we used the bare data of Ref . @ • 
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At tree-level there is a dependence on another factorization scale v s 3> 1/r. The logarithmic dependence 
of V s on v s would cancel the (hard) IR scale dependence of a large but finite mass m (m v s ). Once the 
mass is included, the static energy reads 



E s {r) = 2m(v s ) + V s (r; u a ; v) + 8E s {v) . 



(49) 



The v s dependence of m is the three dimensional relic of the four dimensional pole mass renormalon. In this 
work we only care about the ultrasoft {v) scale dependence. Therefore, we do not need to explicitly consider 
v s , which only appears in the initial matching condition, or, if combined with the mass, is replaced by m 
(up to a constant). On the other hand, note that for AV — V {v) — V s {v) the dependence on v s disappears. 
This is analogous to the ultraviolet renormalon cancellation that takes place in four dimensions, but does 
not mean that AV is renormalon free (in 4D). There is a leftover infrared renormalon, which is reflected by 
the remaining v dependence in three dimensions. 

Note that we have now obtained the exact 0(r 2 ) contribution to the static potential (i.e. the soft contri- 
bution to the static energy). There is nothing left. Moreover, by setting v ~ AT^ potentially large logarithms 
are resummed. 

The above results have been presented in the MS scheme. It is possible to change to a different renormal- 
ization scheme. If we rewrite Eq. (|46|) in terms of AV all the global scheme dependence gets encapsulated 
in the AT^'s and Eq. (|4T|) becomes 



V. 



RG;MS 



{r;v) 



2 2 jr2 



ecu 

13tt 2 



Vy (AV£(r; v) - AV* (r; 1/r)) - %V 2 A r 2 a ( A V| (r; v) - AV» (r; 1/r)) 



2208) C\ 



(19tt 2 + 144) C A n f T F - A8n f T F (4 (tt 2 - 10) C F + ir 2 n f T F ) 



C A A 

x (AV x (r;u) - AV x (r;l/r)) 



384 



(50) 



In a general global renormalization scheme the [AVx(r; 1/r)]™ terms in Eq. (l50l) vanish, when the respective 
matching condition is added, leaving a (globally) scheme independent constant. 

B. 5E%": ultrasoft effects up to NNLL order 

We proceed with computing the static singlet energy. We have already derived the RG improved expression 
for the potential in the previous subsection. Therefore, the computation that is left (if we neglect non- 
perturbative effects) addresses the ultrasoft contribution. The compact expression Eq. (fT2|) includes all 
ultrasoft effects at C(r 2 ) and at any order in the g 2 / AV expansion. At present concrete results are available 
at one, 0(g 2 ), and two loops, 0(g 4 ), given in Eqs. (|13 [ IT4 | . After minimal subtraction they read 



SE[ 



s,MS 



(1-loop) = -aC F Vlr 2 AV^ S I 1 + 7 e - ln(4^ 



2 In 



AV MS 



(51) 



^ s MS (2-loops) = -^a 2 C F AV M sVir 2 ( Ca In 2 



AV MS 



+ C A \n 



AV 



MS 



(■JB ~ ln(4?r)) 



+C A (6 + i ( 7B - ln(47r)) 2 ) - 2n f T F ^ 



(52) 



Setting v ~ AVxiy) resums large logarithms in the potential and minimizes them in the ultrasoft contri- 
butions. Expressed as a double expansion in ar and l/ln(rAV) the result for the static energy reads with 
C(a 2 r 2 ) and NNLL accuracy 



E s (r) = CpaHr 2 ^^) + -C F (7C A - An f T F )a 2 r + V* G ' MS (r;v = AV) + SE^ s MS (v = AV) 



+0{a 4 r z ln u ) 

Ci,aln(r 2 ^ 2 7re^) + ~G f (7Ca - 4n f T F )a 2 
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+C F a i r 2 { -C 2 A \n%rAV MS ) + -C 2 A {2 lE - 1 - In 4) In 2 (r A^ MS ) 



° 2a ( w + 1 ^ + 41r2 2 ~ 7fi(2 + ln256) " 2(11 + ln7r)) 



+n/T F (cU(2 



19tt 2 
48 



^(5-y)j-(-/^) 2 y 



ln(rA^Ms) | + 0(rtV ln u ) , (53) 



where the omitted 0(a 3 r 2 ) terms do not contain logarithms of rAV. Eq. (I5U1) is renormalization scale 
independent up to 0(^^ r ). Furthermore the explicit scheme dependence of V, and 5E^ S and the implicit 
scheme dependence of Eq. (|53|) through the logarithms of AV cancel up to 0(o?r 2 In ). 

Finally we would like to note that the condition v us = AVx{vus), produces a factorization-scale indepen- 
dent scale that is non-pcrturbative in a, 



v us = AVxiyus) = C A aW(l/{C A ad x r)) . 
W(z) is the Lambert function and has the following expansion for large \z\ 

W(z) = ln(z) - ln(]n(*)) - J2(-l) k ln- k (z) £ 



k=0 



(54) 



(55) 



where S 



(-j+k+i) 



is the Stirling number of the first kind 4 . Therefore, v us resums a certain class of logarithms 5 . 



C. Subleading ultrasoft effects: 0{a 3 r 2 ) 

Let us now consider higher order corrections to the result obtained in Eq. (|53"|) . We have already stated 
that, from Eq. (|46l) . the soft contribution to the static energy, i.e. the potential, is known completely at 
0(a 3 r 2 ). This is not the case for the ultrasoft contribution at that order. The general structure of 8E US at 
G(r 2 )is 

SET = C F aV 2 r 2 (AV) 2 £ Y g c n , s ln s . (56) 

n=0 ^ ' s=0 

Interestingly the RG contains important information on the higher order ultrasoft contributions, namely 
on the logarithmic terms. Since the singlet static energy E s is a physical observable and therefore fac- 
torization/renormalization scale independent we find at 0(r 2 ) (neglecting non-perturbative effects for the 
moment) 

, 2 
V ±8E™ = -B{V S ) = -C F aVjr 2 £ B n {C A a) n {AVf- n , (57) 



n=0 



where the B n coefficients can be read off Eq. (|4"5"|) . We obtain the following relation between the coefficients 
at different orders in Ay (note that AV is scale dependent itself): 

B n 8 sfi + (n- 3)c„_ M - (s + l)(c„_ M+ i + c„ jS+ i) = (58) 

with c„ s = 0ifs>n+lorn.<0ors<0 and B n = if n > 2 or n < 0. 



4 see e.g.: http://functions.wolfram.com/IntegerFunctions/StirlingSl 

5 Something similar happens in four dimension when solving the Schrodinger equation with the lowest order Coulomb potential. 
In that case however one sets v B = mCpa(u a ). 
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Ay 




FIG. 2: Schematic tree-level matching procedure for the non-perturbative interaction operator in Eq. f61[) . The 
pNRQCD octet field in the left diagram has been integrated out on the right hand side. 



In particular, Eq. (|58|) allows to fix all logarithmic terms of the three-loop ultrasoft computation: 
SEf™P-< MS (r;v) = C F r 2 a 3 



Ay 



MS 



+ ici(2 + 7B -ln(47r))ln 2 



Ay 



MS 



+ 



C " A ( W - T + \ {lE ~ ln(4 ^ )} + \ {lE ~ ln(47r))2 



+n/2>((2 + ^ )C A + C F (5- £)) - (n f T F f^ 



In 



Ay 



MS 



(59) 



Note that, though we cannot fix the v independent constant C2 t o from RG arguments, it can be computed 
from perturbation theory, but requires a three-loop pNRQCD computation which has not been performed 
yet. The constant term of the static singlet energy at 0(a 3 r 2 ), which is the only missing term in Eq. ([53)1 
to reach N 3 LL precision, is then given by 



y s MS (r;^=-) 



n n 2 2 3 MS 

C F C A r a c 2 



0(a 3 ) 



(60) 



D. Non-perturbative effects: 0{a 4 r 2 /AV) 

At even higher orders in the a/ AV expansion, non-perturbative effects start to contribute. In order to 
study these effects related to loop momenta k ~ a, we integrate out the Ay scale. This means integrating 
out the octet field and ultrasoft gluons. The degrees of freedom left are the singlet field and non-perturbative 
gluons with energy and momentum of order a. The resulting Lagrangian, including the leading order non- 
perturbative effects at 0(r 2 ), reads 

£ np = Trjst (ido - y (r) - 8E?) S j - §^Tr {St( 5 E ■ if S} (61) 

for the case without light fermions (rif = 0). If we were to include light fermions there would also be 
operators of the form 



6£~ ^Tr{tfa 2 qqS} . (62) 

They could generate corrections to the static energy, due to the quark condensate, which are of the same 
parametric order as the purely gluonic ones. In the following we restrict ourselves to the purely gluonic case 
(«/ = 0). 

The coefficient of the non-perturbative operator in Eq. (|6ip is C np = 1 at leading order in the -^V 
expansion. This result is obtained by matching to a pNRQCD tree-level diagram, where two gluons couple 
to the singlet field, as sketched in Fig. [2] 



2 



12 



The interaction with non-perturbative gluons produces a shift of the energy which is proportional to the 
gluon condensate in three dimensions 6 



8E:% = C» O B = 



2tt 



AV b N c (D - l)D 



Ga r^liv 
,.,,11 



(67) 



The leading ultraviolet divergence of the gluon condensate has been calculated in perturbation theory at 
four loops 23 1 . The determination of the finite piece requires lattice simulations [24l l25j| and a computation 
to change from the lattice to dimensional regularization [26[. Taking the result (in the Euclidean) from the 
last reference we have 



b = -(aG^ u G' MU 



) i^Eucliclcan 



-6— C F C 

7T 



43 
12 



157 
768 7 



1 

8^ 



47T 



In- 



2C A g 2 



B G + 0(e 3 ) 



B 



G 



-0.2 ± 0.4(MC) ± 0.4(NSPT) . 



Eq. (|57| is the analog to the Voloshin-Leutwyler correction to the static potential in four dimensions 



■28, 2£ 



which scales like -gy(aG 2 ) 



(aG 2 



In three dimensions we have instead -£y(aG 2 ) 



Renormalizing the bare expression in Eq. (|67|) in the MS scheme yields 



In r \ 



C A C F r 2 a 4 



/43 157 
\~6 ~ 3M n 



157 



In 



-(ln(167r)+ 7B ) 



2B G 



(69) 



To be consistent with our power counting this result has to be combined with the four-loop ultrasoft 
computation. Again the logarithm structure can be determined by demanding scale independence of E s (r), 
which results in the following renormalization group equation for the ultrasoft part 



dv 



-B(V S ) - v—5E np . 
dv 



This equation allows us to fix the coefficients of the four-loop ultrasoft contribution 

C3,4 = , 



C3.3 



-C2,3 



C3,2 = -c 2 ,2 = -- - ^ (lE - ln(47r)) 



c 3,l = C%Y 



r-">' 
L 3,l 



(70) 



(71) 



This is the first term of a Taylor series that can be easily obtained by expanding Eq. (1121) in 1/AV, see also Ref. [St]]: 

oo 

SE n P =J2 G i°i> ( 63 ) 

i=0 

where 

and (H g represents the gluonic Hamiltonian) 

°' = nTd- 1) ( VaC \ E 3 H ? E 3 lvac) ■ (65) 
Using the equation of motion, the gauge fixing and Lorentz covariance we obtain 

O i = Tfg2 v ^...v^v a o...v a ' (66) 

Nl(D-l) 

x<0|Tr([D 8l (0),[...[D /3s (0),G^ op (0)]...][D Q1 (0),[...[D ai (0),G QO ' ) (0)]...])|0>, 

where v is the velocity of the center of mass frame with v 2 = 1 (in the comoving frame v = (1,0)) and the trace is taken in 
the adjoint representation. Generally, however, Eq. 1 163 II does not cover all possible non-perturbative corrections at 0(r 2 ). 
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nUV - 

-3,1 — 



-c-x i represents the coefficient of the ultraviolet logarithm and follows from Eq. §5 



is associated with the infrared logarithm in the ultrasoft contribution, that cancels the renormalization scale 
dependence from the non-perturbative logarithm in Eq. (1691) . 

Thus in order to obtain the complete solution for the singlet static energy up to 0(r 2 a 4 / AV), the miss- 
ing pieces are purely perturbative: C2,o and C3.0, i.e. the non-logarithmic three- and four- loop ultrasoft 
contributions. 

The results given in Eqs. (|53|) . (l59|) and (|69|) have been written as an expansion in 1/ \n(rAV). Alternatively 
one can use 1/ hi(C A ra) as the expansion parameter when setting v us = AVx and using Eq. (|55j) to expand 
AVx in powers of (ln(ln(ra)))™- s / ln s (ra). This expansion might make sense as soon as non-perturbative 
terms like Eq. (f69|) are included, since they can be organized in the same series. Up to 0(r 2 a 2 ) in the 
multipole expansion and N 3 LL order perturbation theory the singlet static energy then reads 



C 2 A C F r 2 a 3 



E s (r) = C F a\n(r 2 v 2 s ire' 1E ) + jC F (7CA-4n f T F )a 2 r 
1 



(73) 



ln J (C A ra) 



+ ln 2 (C A ra) 



ln(CUra) 



1 



~ ln[- ln(CUm)] + h,2 lE - In 4 - 1) 
~ ln[- \n(C A ra)} 2 + ( 7B - In 2) lrx[- ln(C A ra)] + K 



+ - ln[- In(C A ra)] d + -(j E - In 2 + 1) m[- ln(CUra)] 2 + (K + lE - In 2) ln[- ln(C A ra) 



6 



^(^(Y-D-H + SaS)) , fg(7-41n2)(n / T i ,) 2 -C7 F (4- 7 r 2 (|-i^))n / T F 



3C A 



r 12 



X3 

6 



Y-2 209 7T 2 r 715 47 In 2 n 

^-( 7B -ln2 + l)-X(ir + 7£ -l n 2) + — 1 96 24 8 - 



43C(3) 



+ c 2 ,o 



for a general renormalization scheme with cx = ( 7E+ln7r )+ x . 



1 



K = B 2 - cf b + -( 7E + hw)(3 7E - ln(167r)) = 
o 

(2 + ^)n,T F ^(n/IV) 2 + \C F {^ 2 - 10)n f T F 



(74) 



48 

Ca 



a 2 



-( 7£ -ln2) 2 -^ 



11 137T 2 



384 



is a manifestly scheme independent number. Thus the yet unknown ultrasoft constant 02,0 denned by Eq. (|56[) 
is expected to absorb the remaining scheme (X) dependence. 

For rif = we can even go one step further in the l/ln(C A ra) expansion and include the leading non- 
perturbative (Eq. (|69p ) as well as the N 4 LL (ultrasoft four- loop) contribution to the static singlet energy by 
adding 



SE 



C 2 A C F r 2 a 3 



s|N 4 LL 



ln(C A ra) 

+ (K + 7E -ln2 
X 3 1 



\ ln[- HC A ra)] 3 + \{ lE - In 2 + 1) ln[- ln(C A m)] 2 
o 2 



43 
~6~ 



157tt' 
384 



) ln[- \n{C A ra)] - ( 



43 



157tt 2 \ / , , A , 1, 
7B -I- ln(47r) + - 



— - -X 2 ( 7B -ln2 + l) -X X + 7B -ln2+ — - 



6 



43 



G 



384 

1577T 2 

384 



2B G - c 3 ,o 



(75) 



to Eq. (|75|) . The X dependent terms must be subtracted by corresponding terms in the ultrasoft four- loop 
constant 03,0- 

Another convenient choice for the renormalization scale is v = C A a, thus keeping the ultrasoft logarithms, 
but eliminating the logarithms in the non-perturbative finite parts. Eqs. (1731 and (|75|) as parts of an 
expansion up to 0(ln _1 (C A ra)) remain of course unchanged. 
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FIG. 3: Plots of the analytic results for the force F(r) = dE ^ r ^ in ro units with nf =0 in comparison to existing 
lattice data (black dots): 

a) F(r) for N c = 2. The dotted (LO) and dashed (LO+NLO) curves are due to Ref. the results including the 
NNLO order (a 3 r 2 ) in the multipole expansion are new. The labels NNLL and N 3 LL of these curves refer to the 
order in the a/ AV expansion. The most precise result displayed here includes the leading non-perturbative logarithm 
of Eq. $69\) in addition to all corrections up to N 3 LL. We have plotted our results in the MS and the MS scheme in 
order to make the scheme dependence visible and set v = AVx ■ 

b) F(r) for N c = 3. The corresponding curves to a) are shown. 

To further illustrate the relative size of the different orders in the a/ AV expansion, we have displayed the SU(2) results 
for r F(r) starting from the NLO and subsequently adding higher corrections up to N 3 LL including the logarithm of 
Eq. S69V . Plot c) and d) show the situation in the MS and the MS scheme respectively. The LL, NLL, etc. curves 
are understood to be at NNLO precision regarding the multipole expansion. 

Depending on the scheme we have used the values in the first column and the first/third line of tabled for C2,o in the 
N 3 LL plots. 



E. Comparison with lattice data 



The aim of this subsection is, on the one hand, to see how well the short-distance 3D lattice data can 
be reproduced by our theoretical expressions for the static singlet energy and, on the other hand, to try to 
extract numerical values for the ultrasoft three- and four-loop constants C2,o and C3 ; o from fits to this data. 

The authors of Refs. (30l - [3^ | performed quenched (n/ = 0) three dimensional lattice simulations of the 
force F(r) = between two static color sources in the fundamental representation of SU(7V C ) that form 



a color singlet. In Ref. [30( the number of colors is N c — 2 and in Ref. [31( N c = 3. Unlike the static 
energy, the force does not contain the unphysical logarithmic (UV) divergence of the LO static potential, 
which depends on the regularization scheme. It is therefore the preferred quantity to study on the lattice. 
The value of a to be used in our perturbative results for a comparison with the SU(3) data is given by 



7 We don't consider the lattice data for N c = 5 in Ref. [3^1 . because we focus on short distances. 
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a = 2^ [3lj . where r$ — 0.5 /m is the Sommer scale. The corresponding value for a in Ref. [30] was 
obtained from a fit using a wrong ansatz to parametrize the corrections higher than NLO. We repeated the 
fit as described in Ref. [30( using the perturbative result corresponding to Eq. (|53|) in the scheme for which 
d x = 1 (i.e. c x = e 5^ E+ln7r )) and find a = This number is quite similar to the one of Ref. (H (^) 
and indicates that the multipole expansion works quite well for the available SU(2) short distance lattice 
data. 

To illustrate the relative size of the contributions at different orders in the multipole and the ultrasoft 
a/AV expansion, we present two plots of the force F(r) and two plots of rF(r) (normalized to a dimensionless 



number by factors of ro) as a function of r/ro in Fig. [3] We also display the lattice points of Refs. 30] and 31] 
(black dots) for comparison. 

Figs. [3] a) and b) show the results with increasing precision in the multipole expansion for SU(2) and 
SU(3), respectively. At the highest order, i.e. NNLO, there are two curves corresponding to Eq. (|53p in the 
MS and MS scheme with NNLL precision. If we compare the LO, NLO and the NNLO curves with NNLL 
precision, we find a convergent pattern, in particular in the short distance limit. On the other hand the 
a/AV expansion does not work that well, even for the shortest distances that were probed on the lattice. 
Indeed, already at r/r ~ 0.22 for SU(2) and r/r ~ 0.35 for SU(3) we have [s^|] S(7(2) - 0.60 and 

I" jF£ a 1 „,„,, ~ 0.69, respectively. A reasonable convergence of the ultrasoft expansion, say *\j a < 0.25, is 
however only assured for r/ro ^ 0.01 with the given values of a. In the MS scheme it is even worse, since 
generally AVms < AV^g. 

The lack of convergence is highlighted in Figs. [3] c) and d): The LL and NLL contributions are very 
small (maybe anomalously small), much smaller than the NNLL contribution, which, however improves the 
agreement with the lattice data points. The pattern is somewhat similar in the MS and MS scheme despite 
ACL < 2v^s ' ^ e nave a ^ so studied the convergence of Eq. (f73j) . which does not incorporate residual 



AVW ^ AVj 



subleading (N 4 LL) scheme dependences (and is therefore, with the precision of the computation, scheme 
independent). If we consider the equivalent NNLL precision, we find a similar pattern as for the MS and 
MS scheme, though a bit closer to MS. 

We now consider the corrections higher than NNLL in the a/AV expansion. At N 3 LL precision we add 
Eq. (|60p to Eq. (|5"31) in the MS scheme (and accordingly in the MS scheme). In the l/ln(Cyira) expansion, 
the corresponding result is Eq. (|T3|) . At this order the constant C2,o appears in the analytic expressions. It 
is unknown at present. To reach N 4 LL precision we have to add Eq. (|69|) with v — AV to Eq. (|53|) plus 
Eq. (j60|). In the l/h^CUra) expansion, N 4 LL precision corresponds to adding Eq. (|75|) to Eq. (|73|). 8 

From the previous analysis we have seen that the convergence in the a/AV expansion is not good. There- 
fore, we consider any attempt to fit C3.0 as unrealistic and focus on C2,o> for which we only expect an order 
of magnitude estimate from a fit to the available lattice data. 

The strategy we follow is to fit our analytic expressions on the lattice data at the shortest available 
distances in order to obtain 02,0- There is some scheme dependence in this determination, which we use for 
a rough error estimate. On the one hand, we use A): Eq. (|5"3")l plus Eq. (1501) plus the \n{-^-£) term of 

Eq. dnU) for v = AV X in the X =MS and X = MS schemes and, on the other hand B): Eq. ([73]) plus all the 
logarithms of Eq. (|75[) . Ideally most of the scheme dependence should originate from the X dependent terms 
in Eq. (|73[) . We also vary the number of lattice data points in the fit to see the dependence of the results. The 
numerical value should not depend much on the number of data points if we are at short enough distances. 
Finally, we also consider determinations for N c = 2 and N c = 3, as we expect C2,o to be independent of N c 
(for n f = 0). 

TableUlists the numbers for €2,0 resulting from the fits to the data points at the smallest available distances 
r in the renormalization schemes MS and MS. Albeit with large uncertainties, we see that the numbers are 
in the same ballpark. We should moreover not equally weight all the determinations. We believe for instance 
that the MS fit (A) is less reliable, since AVms < AV^jg. We also expect that increasing the number of 
included data points (with larger r) decreases the fit quality. In any case, we should not forget about the 
poor convergence, which makes a precise determination of C2,o impossible. On the other side we find it 
reassuring and nontrivial that the SU(2) and SU(3) determinations, just as the (A) and (B) MS results, are 



In the following we will only take into account the logarithmic terms in Eqs. J69H and H75I I. since the constant term includes 
C3,0. 
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gauge group: 


SU(2) 


SU(2) 


SU(2) 


SU(3) 


# of fit points: 


1 


3 


5 


1 


(A): MS fit 


-0.04 


-0.42 


-0.82 


+0.41 


(A): MS fit 


+0.92 


+0.18 


-0.67 


+0.85 


(A): MS from MS fit 


+2.64 


+2.26 


+1.87 


+3.09 


(A): MS from MS fit 


-1.76 


-2.51 


-3.35 


+0.18 


(B): MS fit 


-0.74 


-1.22 


-1.75 


-0.43 


(B): MS fit 


+1.95 


+ 1.47 


+0.93 


+2.25 



TABLE I: Results for C2,o from fitting the force derived from (A); Eq. S53p plus Eq. f60\) plus the logarithm of 
Eq. fgPp, and (B); Eq. (75p plus all the logarithms of Eq. (75$ , to lattice data for N c = 2 and N c = 3 in the X = MS 
and X = MS scheme. The (first/first three/first five) lattice points with smallest distances r have been used for the 
fit. To obtain the values in the third and the fourth line we have transformed the MS fit values to the MS scheme, and 
vice versa, by adding /subtracting the difference coming from the scheme (X ) dependent terms in Eq. (7ff| ). Note that 
for option (B) such a scheme transformation reproduces the same results as from the direct fits. There is therefore 
no uncertainty from the scheme dependence in this case. 



quite close. We have also performed fits not including the non-perturbative logarithms. Their effect is small 
and within the uncertainties, which, from the above discussion, we roughly estimate to be as follows: 

cff~ -0.04+2. (76) 

We took the central value from fitting only to the first lattice point at r/r Q ~ 0.22, where we expect the 
smallest theoretical error from the truncation of the perturbation series. 

We use this value for the MS N 3 LL and N 3 LL+non-perturbative-logarithm curves shown in Figs. [3] 
In other words we have fitted the N 3 LL+non-perturbative-logarithm curve, which includes all mentioned 
corrections in the MS scheme, to the first data point for N c = 2. That's why this point lies exactly on the 
line, whereas in Fig. [3jb), where N c = 3 and the same value for c^q nas been used, it does not. As we 
have already stated the convergence is not good. This is also seen in the difference between the NNLL and 
N 3 LL curves, although the difference between the N 3 LL and the N 4 LL curves due to the non-perturbative 
logarithms is small at short distances. The pattern in the MS scheme is similar even though the convergence 
is worse. Working with Eq. (|73|) plus all the logarithms of Eq. (|75l) produces curves similar to those obtained 
in the MS scheme. 

From the above analysis, we conclude that the lattice simulations so far have produced data at short 
enough distances to test the multipole (or) expansion (cf. Figs. [3] a) and b)). A Precise test of the a/AV 
expansion is however not possible with the existing data. 



IV. PNRQCD(D=2) 

It is also interesting to investigate the static quark-antiquark potential in two space-time dimensions. In 
the following we set rif = 0. 

For D = 2 the exact result for the singlet static potential is known. It can be derived from its standard 
(non-perturbative) definition via a Wilson-loop: 

E s (r) := lim ± In (-^TrP exp ( - ig I A^dx*)) , (77) 

where the contour T of the integral in Eq. ([77]) is a rectangle with spatial extension r and extension T along 
the time coordinate axis in Minkowski space. Choosing axial (A\ = 0) gauge in exactly two dimensions, the 
only non- vanishing components of the gluon field-strength tensor are i*oi = — -Fio = —9iAq. Hence, (Aq) 
gluons do neither interact among themselves, nor propagate in time, since no time derivative is left in the 
Lagrangian. Therefore only planar "ladder" diagrams (with "potential" gluons) contribute to the Wilson- 
loop. In such a quasi- Abelian situation, one can forget about the path-ordering operator V in Eq. (|77|) and 
treat the SU(JV C ) generators T a in the exponential as if they were proportional to the unit matrix. Finally 
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we arrive at the result for the corresponding Abelian Wilson- loop multiplied by jj-Tr T a T a — C F . This 
result is identical to the perturbative tree-level contribution 

E s (r) = 2nC F ar (78) 

and therefore proving exponentiation. 

The computation in pNRQCD is organized differently, because of the factorization of the different scales 
that contribute to the potential. For exactly D = 2 pNRQCD recovers Eq. (|75)) . since the soft contribution 
is precisely Eq. ([75)1 and there are neither ultrasoft nor non-perturbative gluons. For D > 2 one would 
expect the result to converge to Eq. (|78|) . when D — > 2, provided the limit is smooth. The soft two-loop 
calculation [2j yields however 

V s (r) = 2nC F ar - 2irC A C F a 2 r 3 + ^(33 + 2ir 2 )C 2 A C F a 3 r 5 (79) 

for D —¥ 2 and rij — 0. From the effective theory viewpoint, we expect that the ultrasoft contributions to the 
static energy, calculated in pNRQCD for D — > 2, cancel the a 2 and a 3 terms in Eq. (|79|) . Though, checking 
this is a delicate issue, since as mentioned in subsection lll Cl the proper ultrasoft expansion parameter is now 
A¥- ~ gr, i.e. non-perturbative. Thus the intrinsically non-perturbative ultrasoft (AV/ g) expansion and the 

perturbative multipole (gr) expansion mix. 

The non-perturbative nature of the ultrasoft expansion does not allow a calculation in terms of pNRQCD 
loop diagrams for I) / 2. Expressions like Eq. (fT2j) should rather be expanded for small AV and the full 
non-perturbative result for the chromoelectric correlator should be used. Of course the latter is not known 
in dimensional regularization. 

Despite that we have observed interesting cancellations. The ultrasoft one-loop contribution at the lowest 
non- vanishing order in the multipole expansion, Eq. (|13p . gives 

5E us (l - loop) = 2irC A C F a 2 r 3 (80) 

in the limit D — > 2 and exactly cancels the soft a 2 contribution in Eq. (|T9"1) . 

At this point we are unfortunately not able to check the (anticipated) cancellations at higher orders in a, 
as they may require higher order terms in the multipole expansion. We are also not sure about the relevance 
of the above result, as it follows from a perturbative ultrasoft computation. Nevertheless the observed 
cancellation could be interpreted as another check of the applied effective field theory approach, which uses 
perturbative methods to derive meaningful results for the singlet static energy in arbitrary dimensions. 



V. CONCLUSIONS 



The main achievement of this paper is a precise determination of the static energy of a color singlet 
quark- antiquark system in three space-time dimensions for ar < 1. We have complemented the strictly 
perturbative (soft) two-loop QCD result Q by an effective theory calculation in pNRQCD to account for 
effects at the ultrasoft scale (AV — V — V s ), which are not accessible to pQCD. 

We observe a perfect cancellation of the IR singularities in the soft pQCD result at 0(a 3 r 2 ) and the UV 
divergences in the one- and two-loop pNRQCD corrections, i.e. at NLL and NNLL order in the ultrasoft 
a/ AV expansion. This represents an independent and nontrivial check of the soft pQCD and the ultrasoft 
pNRQCD contributions, because the singlet static energy is a physical observable and must be IR finite. We 
also observe a cancellation of logarithmic divergences that can be associated with the leading renormalons 
of the potential and the quark mass in four dimensions. 

Demanding complete IR finitencss we have determined indirectly the UV divergences at ultrasoft three- 
loop level, without having calculated the actual diagrams. In fact the only missing piece in our N 3 LL 
prediction for the singlet static energy up to 0(r 2 ), Eq. (f73|) . is the ultrasoft three-loop constant ca,o- We 
also provide expressions for the N 4 LL correction and the known leading non-perturbative contribution from 
the gluon condensate 23-26|, which is parametrically of the same order. They are complete up to the 
ultrasoft four- loop constant c^^q. Thus the precision is currently limited by missing perturbative three- and 
four-loop computations. 

From a fit to SU(2) lattice data [3(| we were able to extract a rough value for C2,o. The uncertainty is 
however rather large, because data points for distances r sufficiently small to ensure reasonable convergence 
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of the perturbative a/AT^ series are currently not available. The ar (multipole) expansion in contrast can 
be tested with the present precision of the lattice simulations. 

With lattice data at smaller distances and analytic results for c^n and 03,0, an interesting cross-check of 
the three-dimensional lattice simulation for the gluon condensate [24], HE] in Eq. (1551 would be feasible. 

We have also studied the two dimensional case, for which the exact result is known and equivalent to the 
tree-level result. Therefore, strong cancellations have to occur among the contributions from the different 
scales: soft, ultrasoft and g in the D — > 2 limit. We have found a cancellation between the one-loop soft and 
the one- loop ultrasoft results in dimensional regularization for D — >■ 2. The relevance of this cancellation 
is however not clear to us, because it follows from a perturbative ultrasoft computation. Nevertheless, we 
think that it could be worthwhile to explore the D — > 2 limit in more detail, and see whether it could provide 
nontrivial information on the structure of the perturbation series for general dimensions. We leave this issue 
open in this paper. 
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